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Abstract 

The mean-field theory with the use of the slave-boson functional method is 
generalized to take account of the spin- and orbital-ordered state in the dou- 
bly degenerate Hubbard model. Some numerical calculations are presented 
of the antiferromagnetic orbital-ordered state in the half-filled simple-cubic 
model. The orbital order in the present theory is much reduced compared 
with that in the Hartree-Fock approximation because of the large orbital fluc- 
tuations. From a comparison of the ground-state energy, the antiferromag- 
netic orbital state is shown to be unstable against the antiferromagnetic spin 
state, although the situation becomes reversed when the exchange interaction 
is negative. 



*To be submitted to a special issue of "Foundation of Physics" celebrating the 75-th birthday of 
Martin Gutzwiller 

t E-mail : hasegawa@u-gakugei .ac.jp 

1 



I. INTRODUCTION 



In his seminal paper, Gutzwiller |]] proposed a variational approach, employing the pro- 
jected wave function to take into account the effect of correlated electrons. Since the exact 
evaluation of the ground-state energy is difficult, he introduced an additional approximation 
which is now called the Gutzwiller approximation (GA). A validity of the GA was studied 
by Monte-Carlo method for finite-size clusters [[J ||. It is realized that the GA becomes a 
better approximation for higher- dimensional systems, and it is exact in the limit of infinite 
dimension. The GA has been widely adopted in many area including metal-insulator tran- 
sition (MIT) HI, magnetism and high-T c superconductors (for a review of the GA of SHM 
see Ref. [|). 

Kotliar and Ruckenstein developed the slave-boson mean-field theory by using the 
saddle-point approximation to the slave-boson functional-integral method. They successfully 
derived the antiferromagnetic solution besides the paramagnetic solution of the GA. It is 
shown that the GA and the slave-boson mean-field theory is equivalent. 

Most of the theoretical studies based on the GA have been made for the single-band 
Hubbard model (SHM) for its simplicity. Actual systems, however, inevitably have the 
orbital degeneracy. It is necessary to investigate the role of the orbital degeneracy and the 
effect of Hund-rule coupling due to the exchange interaction for a better understanding on 
strongly correlated systems. 

The first attempt to include the orbital degeneracy within the GA was made by Chao 
and Gutzwiller ||. In the last few years the Hubbard model with orbital degeneracy has 
been extensively studied by using not only the GA-type theories [0- |l4j] but also the dy- 
namical mean-field approximation ||15||, []16|], and Monte-Carlo simulations W7\- |]19||. The 



original GA proposed by Chao and Gutzwiller || was reformulated in Refs. 0- [10]. Lu 
|7]] obtained the analytical expression of the critical interaction for the MIT. Okabe [§] pro- 
posed a sophisticated method in calculating te band-narrowing factor which is a difficult 
part in applying the GA to the degenerate band model. The first-order transition with an 
introduction of the exchange interaction was discussed in Ref. |J. Biinemann, Weber and 
Gebhard |T(| made a detailed study by using the generalized interactions. 

The present author proposed a slave-boson mean-field theory for the degenerate 
Hubbard model, by employing the slave-boson method proposed by Dorin and Schlottman 
PHI for the Anderson lattice model. Fresard and Kotliar [TJ] adopted an alternative slave- 
boson functional integral method. The slave-boson mean-field theory of Hasegawa ]ll]] and 
Fresard and Kotliar [|nj] is the simple generalization of the Kotliar-Ruckenstein theory PJ 
for the SHM to that for the degenerate Hubbard model, and it is again equivalent to the 
GA 0- ]l(|. The MIT in the paramagnetic state in the doubly degenerate Hubbard model 
(DHM) is discussed ]13| |T^] . Subsequently the slave-boson theory was extended and applied 
to the antiferromagnetic (spin) state of the half-filled DHM [12| . 

One of the advantages of the slave-boson mean-filed theory to the GA is that it has the 
wider applicability than the GA. By using the Green-function method, for example, the 
slave-boson mean-field theory can easily deal with the antiferromagnetic state in the SHM 



21 1, [22] and in the doubly degenerate Hubbard model (DHM) [0. This suggests that 



we can discuss a system with more complicated orderings with the slave-boson mean-field 
theory. It is the purpose of the present paper to generalize our theory to include the orbital 
ordering besides spin ordering. 
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The paper is organized as follows: In the next Sec. II, we present a basic formulation of 
our slave-boson saddle-point approximation to deal with the spin- and orbital-ordered states. 
Since calculations for the paramagnetic and antiferromagnetic states have been reported in 
our previous papers |IT|| - H13] , we present, in Sec. Ill, some calculated results for the orbital- 



ordered state in the half-filled DHM. Section IV is devoted to conclusion and supplementary- 
discussion. 



II. FORMULATION 

2.1 Basic Equations 

We adopt the Hubbard model with the doubly orbital degeneracy whose Hamiltonian is 
given by 

# = EEE ^ r '4.c JT v + ~ E E ctfcLcwi^cw 

a ij tt' i (t,ct)^(t' ,a') 

_ EE E( a/l i + T 9i) claCiT'a, (1) 

air 

where Ci Ta is an annihilation operator of an electron with an orbital index r(= ±1) and spin 
a (= ±1) on the lattice site i. The first term expresses electron hopping, which is assumed 
to be allowed only between the same subband: t\J = Uj5 TT >, for a simplicity. The on-site 
interaction, U™, , in the second term is assumed to be given by 

U% = U = U for r = r', a ^ a', (2) 
= U 1 = U-2J for r ^ r> ^ a', (3) 
= U 2 = U - 3J for r ^ r', a = a' } (4) 

where U and J are Coulomb and exchange interactions, respectively. In the third term of 
Eq.(l), hi and ^ are the magnetic and crystal fields, respectively, at the i site. 

In order to discuss the DHM with the slave-boson functional integral method, we in- 
troduce, for a given site i, sixteen slave bosons which are classified into the following five 
categories: 

(1) ej for the empty state; 

(2) p iT for the singly occupied state with a cr-spin electron in the r band; 

(3) dior for the doubly occupied state with a pair of up- and down-spin electrons in the r 
band; dn T for that with an up-spin electron in the r band and a down-spin electron in the 
— r band; d^a for that with a-spin electrons in the two subbands; 

(4) t ira for the triply occupied state with a pair of up- and down-spin electrons in the r 
band and an extra a-spin electron in the — r band; and 

(5) fi for the fully occupied state. 

These slave bosons obey the constraint given by 

e i + E Pircr + E(^°r + d ilr ) + ^ d i2a + E E **nr + fi = 1, (5) 

IT(T T <J <J T 

and the equivalence between fermion and boson operators is given by 
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^ira^iTa Pira d{0 T -\- dn aT -\- A%2cr ~\~ tvra tir—a ^i—ra fi 



(6) 



By incorporating the conditions given by Eqs. (5) and (6) with the Lagrange multipliers, 
and \fj) a , we get the partition function given as a functional integral over the coherent 



state of Fermi and Bose fields 11 



Z = jDejDp J Dd J Dd! J Dd 2 J Dt jDf Jd\ {1) JdX^ exp - J (L f {t) + L b {t)) 



(7) 



with 







L f( f ) = EEE Z ira Uj Zjra 4ra C jra + EEE(^ + A Sr ~ * K ~ T 9i) c\ T(7 C iTU . (8) 



i.l 



l <J T 



Of 



L b (t) = £ 



J2(U d l0T + U x d ilT ) + U 2 J2 d i2c + (U + U x + U 2 ) ((£ E Ut*) + 2 fA 

t a \ a r / 

e i + £ V%tg + J2(d i0r + d ilT ) + £ d i2a + £ £ ^™ + fi ~ 1 







0-h 



^ ] ^ ] ^ ] ^ira \Pircr <^iOr d{\ aT ~\~ d{ 2a -\- tj T(7 ~\~ tj T _ a ~\~ ti^ Tcr -\- /j 



(9) 



i a t 



Here Dp = njlLj-Il,- d p iTU et al, f3 is the inverse temperature, Lf and L b denote the terms 
relevant to the Fermi and Bose fields, respectively, and z ira is given by 

Ziro (^ircr) ^ i^J^iPir a \/Pir —crdiQ T ~\~ "s/pi—r—adiXar \fpi—r a di 2a 

+a/ dio -r ti_ Tcr + y/ dn - aT ti Tcr + yj di 2 _ CT ti T _ CT + \J di_ T — (j fi) (1 — Tii Tcr ) . (10) 

where n ira = < c\ Tcr Ci Ta > and the bracket denotes the expectation value. 

In discussing the spin- and/or orbital-orderded states, we make the following change of 
variables for n; TrT : 



^ = ££ 



11; 



(11) 



M, = EE^™ (12) 

a t 

°i = EE rn ^' (1 3 ) 

a r 

p * = EE arr ^ CT , (14) 

a t 

where Ni, Mi, Oi and Pi denote operators relevant to the number of electrons, the magnetic 
moment, the orbital order and the orbital spin polarization, respectively, at a given % site. 

(2") 

Similarly, we make the change of variables for \ iTa as follows: 
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^ = (1/4) EE«. ( 15 ) 

<T T 

6 = (-1/4) EE^ A £, ( 16 ) 

<t r 

^ = (-1/4) EE (17) 

a t 

* = (-m EE^€> (is) 

it r 

where i/j, £^ 0j and stand for conjugate fields of N iy Mi, Oi and Pi, respectively. Then 

ssed by 

n ira = (1/4) (N t + a Mi + TO t + or Pi) (19) 



Hirer and \ f T \ are expressed by 



All = iyi- o ii - t (\>i - or rji) (20) 

Substituting Eqs.(19) and (20) to Eqs.(7)-(9) with some manipulations, we get the par- 
tition function within the static approximation, given by 

Z = JDN J DM J DO J DP J Dv J Di J Dcj) J Dn 

x J Dd j Ddi j Dd 2 j 'iX j Df exp(-/3$), (21) 

with 

$ = $ + $i + $ 2 , (22) 



$o = E 



EWiOr + U.dar) + U 2 E d i2u + (C/ + C/i + U 2 ) f(E E *ira) + 2 /») 

. r a \ a t ) _ 



(23) 



$i = E [(^ - ^ + + + ^ > ( 2 4) 



$ 2 = y cfe f( £ ) (— 1/tt) Im Tr lnG(e). (25) 

Here /(e) is the Fermi-distribution function, e F the Fermi energy and G(e) is the one-particle 
Green function given by 

G(s) = (s-H cS y 1 , (26) 
where the effective slave-boson Hamiltonian, H c s, is given by 
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H oS = E E E <l l Tatijct T(T c jTa + E E N - a ( h i + &) - T (9i + <t>i) - rii] c\ Ta c iTa . 

ij a r % a t 

(27) 

and the band-narrowing factor, q l J a , is given by 

Qra = Ziro 2jV(7) (28) 

with 

Zircr Ci T<T / Di T<T) (29) 



C iTa = lQ(^/eiP iTa + ^Pi T -adiO T + yjpi -t —o-dii aT + \J Pi—T a di2 a 



+ V diO -r ti- T a + \Jdi\ - ar ti Ta + \J di2 ~cr U T -a + \J d{- T _ CT /j), (30) 

Di Tcr = (Ni + aM t + t0 1 + t aP,) 1 ' 2 (A-Ni- aM, - rO z - r aP,) 1 ' 2 . (31) 

d = 1 - Ni + ^(d i0T + d ilT ) + d i2a + 2 E E Ura + 3/, (32) 
T a a t 

p iTU = (1/4) (AT* + <j Mi + r Oi + arPi) - (d i0r + d n aT + d i2a + E + *»-t<t + /i) (33) 

Since the effective transfer integral in Eq.(27) is expressed as a product form: 
ZiTatijZj Ta , we can express G(e) and $2 in terms of the locator defined by 

Xira(s) = [e-i/i + a(hi + &) + r(g + (pi) + arrji] /r iTa , (34) 

where 

r iTU = (Zircr) j (35) 

explicit expressions for $ 2 being given shortly [Eqs. (61), (62)]. 

By evaluating the saddle-point value of $, we get the slave-boson mean- field free energy 
given by 



J2(Uod l0 r + U^t) + U 2 J2 d i2 * + (Uq + C/l + £/ 2 ) ( (E E **nr) + 2 /j) 

r (T \ a t / . 

+ E - "i) N i + + <t>i°i + vtPi] 

i 

+ J de f(e) (-1/tt) Im Tr lnG(e). (36) 

It is noted that for a given set of parameters of U, J, T and N, the total number of electrons 
(per site) given by 
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N= (l/JVb)EJVi = (l/iV )EEE^-. (37) 

i i a t 

where N is the number of lattice sites. The optinum values of iVj, Mj, 0«, P«, i^, £j, fa, 
?7i, djor, <^iir, ^i2o-, U Ta , fi and £f in Eq.(36) are determined by the self-consistent equations 
given by Eqs.(ll)-(14) and the following equations: 

£F-"i + EE RirAdTirJdNi) = 0, (38) 

<T T 

6 + EE RirAdTirJdMi) = 0, (39) 
<t r 

0i + E E Rira(dr iT<T /dOi) = 0, (40) 

<T T 

^ + E E Rirvidr^/dPi) = 0, (41) 

ct r 

f/o + EE Rir>v{dr iT , a /dd iGT ) = 0, (42) 

t/i+EE RiM9nr'a/dd ilT ) = 0, (43) 

f/ 2 + EE RiT*>{dr iTa ,/dd i2a ) = 0, (44) 

C/ + C/i + C/ 2 + E E ^ {dr iT , a ,/dt iTa ) = 0, (45) 

a' r' 

2(C/ + U! + U 2 ) + E E RiraidTirJdU) = 0, (46) 



with 



^ = d$ 2 /dr iTtr . (47) 



n ircr = J de f(e) (—1/tt) Im < ira \ G(e) \ ira >, (48) 

where r ira is given by Eq.(35), and the bra < ira | and ket | ira > denote the cx-spin state 
in the r subband at the % site. 

2.2 Spin- and Orbital- Ordered States 
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The expression for $2 given by Eq.(25) (and then Ri Ta in Eq.(47)) depends on G(e) 
given by Eqs.(26) and (27) which is specified by the electronic structure and by the spin- 
and orbital-ordered states to be investigated. In order to discuss the spin- and/or orbital- 
ordered states in the paramagnetic, ferromagnetic or antiferromagnetic state on the same 
footing, we divide the crystal into two sublattices, A and B. We assume that for the AF 
wave vector Q, the relation: Eu+q = —Sk holds where Sk is the Fourier transform of the 
transfer integral, tij. 

The locator given by Eq.(34) at the i (j) site belonging to the A (B) lattice, is given by 



Aj r(7 



(e) = X At(T (e) = [e -v A + a(h A + 60 + r(g A + <p A ) + <jtt) A \ /r Ari 
X jTa (e) = X BTa (e) = [e-u B + a(h B + f B ) + r(g B + <f) B ) + arr] B ] jr Bl 



(ieA), (49) 
(jeB). (50) 



When the spin and orbital orderings in the A and B sublattices have the simple symmetry 



relation such as v A = v B = v, £ A = ± £ B = £ and r] A = ± r] B = rj with h A 
b B = 0, Eqs.(49) and (50) become 



X At<t {e) 
X Bt Je) 



e - v + cr£ + T(f) + crrr])/r Ar 
e - v + cr£ + r0 + o-rr])/r BT 
e — v — cr£ + T(j) — arr])/r BT 
e — v + <j£ - rcj) — o-Ti])/r BT 



= (e — v — <t£ — tc, 
with r Aru and r Bro . given by 

= (^tct) ) 



(TTT])/r B7 



= r Aj 
= r Bj 



for F-f state, 
for AF-f state, 
for F-af state, 
for AF-af state, 



(ieA) 



= f 'Atcf 
= r Ar-cr 
= r A-rtj 



for F-f state, 
for AF-f state, 
for F-af state, 
for AF-af state. 



U G B) 



= h B = g A = 

(51) 
(52) 
(53) 
(54) 
(55) 



(56) 
(57) 
(58) 
(59) 
(60) 



In Eqs.(51)-(60) "AF-f state", for example, expresses the state in which spin and orbital 
orders are antiferromagnetic (AF) and ferromagnetic (f), respectively [when spin (orbital) 
order is vanishing, it is the paramagnetic state and is referred to as P (p) state]. 
After a simple calculation, $ 2 in Eq.(25) becomes 



$2= fdef(e) (l/7r)Im £EE In (q. 



2 

rcr 



X ATa (e)X Bro .(e) — e 



(61) 



with the band-narrowing factor, q Ta , given by 

<?rcr = \Jf' AT(J r Bl 



(62) 



where r Arc and r Bro . are given by Eqs.(56)-(60) depending on a kind of the spin and orbital 
orderings. The expressions for R iTa and n iTa given by Eqs. (47) and (48) become 



R lTa = J def{e) (-1/tt) lm[(Q Ta /r iTa ) F (i2 T 



(63) 
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rw = de f(e) pi T a(e), (64) 



where pi Tcr is the local density of states of electrons with spin a in the band r on the site i, 
given by 

pira(e) = Im [K ATa (e)/r A ra] (i e A), 

= (-1/tt) Im [K Br(T ( £ )/r Br(7 ] G 5), (65) 

with 

K Ara (e) = [X BTa (e)/X Ara (e)} 1/2 F (Q Ta ), (66) 
K BTa (e) = [X ATC {e)/X BTC {e)\ 112 F (O r(T ), (67) 
fi T<T (e) = [X Ara ( £ )X BT(T ( £ )] 1/2 , (68) 

F (e) = Jde' Po (e')/(e-e'), (69) 

being the unperturbed density of states. The optimum values of fields of u, £, ip and 
rj, the order parameters of M, O and P and the occupancies, the average of the bosons of 
d 0r et ai, are self-consistently determined by Eqs.(ll)-(14) with Eqs.(38)-(48). 

The expressions for the free energy and the self-consistent equations given by Eqs.(36) 
and (38)- (48) without the orbital ordering (O = P = and = rj = 0) reduce to those 
of the GA for the P-p state @, || and to that for the AF-p state Ul2j . Thus the present 
slave-boson mean-field theory becomes the generalization to the cases including both the 
spin- and orbital ordered states in the DHM. 

2.3 Comparison with the Hartree-Fock Theory 

It is worth to compare the expression for the free energy in the present theory given by 
Eq.(36) with that in the Hartree-Fock approximation given by 

F HF = (1/2) [fa ~ Vi)Ni + (iMi + frOi + rtiPi] (70) 

i 

+ Jde f(e) (-1/tt) Im Tr \nG HF (e). 

with 

G BS {e) = {e-H m )- 1 (71) 

HrF = EEE tijclraCjra + N _ + &) _ T (di + <Pi) ~ °T Tfr] 4 Ta C iT(7 , (72) 

ij a t 

u i: 4>i and rji being given by 



l <T T 
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i/< = (l/4)(Uo + Z7i + Z7 2 )JVi, 



(73) 



&=(l/4:)(U + U 1 -U i )M i , (74) 
fa = (l/A)(-U + U 1 + U 2 )O l , (75) 

7i i = (l/4)(y -U 1 + U 2 )P i . (76) 

Here iVj, Mj, Oi, Pi and are determined self-consistently by Eqs.(ll)-(14) with (73) -(76). 
It should be noted that £j, fa and rji in the Hartree-Fock approximation are given by 
Eqs.(73)-(76) while those in the present theory are variationally determined by Eqs.(38)- 
(41). The free energy of our theory is expected to be generally lower that the Hartree-Fock 
free energy, as numerically shown for the P-p, AF-p and P-af cases. 

III. CALCULATED RESULTS 

The formalism presented in the preceding section can be applied to various types of 
spin- and/or orbital-ordered states. Since calculated results for P-p and AF-p states have 



been published in our previous papers [11]- [13], we here present some calculations only for 



orbitally ordered (P-af) state in order to demonstrate the feasibility of our theory. 

We perform numerical calculations for the simple-cubic model with nearest-neighbor 
hoppings. Input parameters for our calculations are the non-interacting density of states, 
Po(e), the Coulomb and exchange interactions, U, J, and the total number of electrons, 
N, which is two for the half-filled case under consideration. We employ the approximate, 
analytic expression for po(e) of the simple-cubic density of states |T2[. The energy and the 



interactions are hereafter measured in units of a half of its total band width. The calculated 
ground-state energy without interactions (U = J = 0) is e = —0.3349, which is in good 
agreement with the exact value of —0.3341 [RJ. Details of the calculation method is explained 



in Ref. 12 



A. J = Case 

We firstly discuss the case of vanishing exchange interaction (J = 0), for which the P-af 
state is equivalent to the AF-p state discussed in Ref. fli~2"] . The U dependences of the orbital 
order O, band-narrowing factor q, and the energy E of the P-af state are given by those 
of the sublattice magnetization M, q, and E of the AF-p state shown in Figs.l, 3 and 5, 



respectively, of Ref. [12]. The U-dependence of the occupancies of the P-af state is obtainable 



from that of the AF-p state shown in Fig. 4 of Ref. |1"2]]; do a = dob = do, d\ a = dxb = d\ 



4t = <4 T = d h d al = d bl = d b t aT = t feT = t h and t al = t H = t l in Fig.4 of Ref. JT2 
should be read with the following changes: d oa <-> d 2 -[, g? & <-> f%, t a j. — > t&j, — > t&j, and 
t^} tai, ^ ne others being the same for the AF-p and P-af states (a and b denote the two 
orbital bands). 
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B. Finite J Case 



Next we introduce the exchange interaction J into our calculation. In the case of finite 
J, the equivalence between the AF-p and P-af states mentioned above for J = is no more 
preserved. Solid curves in Fig. 1 show the J dependence of the orbital order O for [7=0.5, 
1.0 and 1.5. For a comparison, results of the HFA are also shown by dashed curves. With 
increasing the value of J, the orbital order O in the P-af state decreases. Because of the 
large orbital fluctuations, the orbital order in the GA is much smaller than that in the HFA, 
and it almost vanishes for J > 0.05 where the P-p state is realized. This is in contrast with 
the sublattice magnetization M in the AF-p state, which increase as increasing J ||12|| . This 
is easily understood because the effective fields of O and M are given by (U — 5J)/4 and 
(U + J)/4, respectively, in the HFA (Eqs.(2)-(4),(73),(74)). 

The densities of states of the P-af state for J = 0.02 and 0.0 with U = 1.0 are shown 
in Figs. 2(a) and 2(b), respectively, where p a ^ = p a [ = p a and = p^ = pb- It is noted 
that Fig. 2(b) for J = expresses also the density of states of the AF-p state if we read 
Pa - > Pt( = Pal = Pbi) an d Pb ~ ► Pi(— Pal = Pb[) [0 • When the J value is increased from 0.0 
to 0.02, the orbital order decreased and the polarization in the densities of states decreases 
as shown by Figs. 2(a) and 2(b). 

Figure 3 shows the J dependence of the band-narrowing factor q for [7=0.5, 1.0 and 1.5. 
With increasing J, q increases slightly for [7 = 0.5 while it decreases for [7= 1.0 and 1.5. 

The J dependence of the occupancies is shown in Fig. 4, where from the symmetry, we 
get c?2t = <^2| = <^2, t a ] = t a i = t a and tb] = hi = h- When J is increased, d Q b, di a , du, 
t a and / increase whereas d 0a and dob decrease. At J > 0.05, the system becomes the P-p 
state where do a = dob and t a = tb (see Fig. 1). 

In order to investigate the relative stability between the AF-p and P-af states, we calcu- 
late their total energies E, which are shown in Fig. 5. The bold (thin) solid curve expresses 
E of the AF-p state in the GA (HFA), and the bold (thin) dashed curve E of the P-af state 
in the GA (HFA). For J = 0, the AF-p and P-af states have the same total energy. We note 
that for all the J values investigated, E in the GA is always lower than that in the HFA for 
the Af-P and p-AF states, as well as for the P-p state [[T"2]j . For finite (positive) J, the 
energy of the AF-p state is lower than that of the P-af state both in the GA and HFA. 



IV. CONCLUSION AND DISCUSSION 

To summarize, we have developed the theory of the spin- and orbital-ordered states in 
the DHM, generalizing our slave-boson mean- field theory [[II]]. In order to demonstrate the 
feasibility of our theory, we have presented numerical calculations of the antiferromagnetic 
orbital (P-af) state, showing the J dependences of the orbital order O, the band narrowing 
factor q and the occupancies. From the calculation of the total energy E, the antiferromag- 
netic orbital (P-af) state is shown to be unstable against the antiferromagnetic spin (AF-p) 
state except for J = for which the both states are energetically degenerate. 

In order to more investigate the role of the exchange interaction on orbital order in the 
DHM, we have made numerical calculations also for the negative J, although J is conven- 
tionally taken to be positive. We notice in Fig.l that for the negative exchange interaction, 
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the orbital order O is considerably increased: we cannot obtain the solution at J < 0.048 
(J < 0.014) for U = 1.0 (U = 1.5). This increase in O is induced by a significant increase 



in do a , as shown in Fig. 4. As was discussed in Ref. [0, the negative exchange interaction 
is not favorable for the sublattice magnetization M in the AF-p state. Then the increase in 
O and the decrease in M for the negative J lead to the P-af state to be more stable than 
the AF-p state, as the energy calculation shows in Fig. 5. 

We previously showed from a comparison between the calculations using the GA EB 



and the variational Monte-Cairo method || that the GA is a fairly good approximation for 
the one-, two- and three-dimensional simple SHM (the GA is exact for infinite dimensional 
SHM) . We expect that this would hold also for the DHM although no variational calculations 
for the DHM have been reported yet. We can apply our theory to states in which both spin 
and orbital orderings coexit such as the AF-f and F-f states. A study of the N-U phase 
diagram (N: electron number, U: interaction) of the simple-cubic DHM, as was made within 
the HFA ||24j| , will be interesting. It would be also promising to investigate the temperature- 



interaction phase diagram of the DHM, by generalizing our approach [^] [^5| in which the 



effects of electron correlation and thermal spin fluctuations are properly taken into account. 
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FIGURES 



FIG. 1. The J dependence of the orbital order O in the P-af state calculated by the GA (solid 
curves) and HFA (dashed curves) 

FIG. 2. The local densities of state in the P-af state for (a) J = 0.02 and (b) J = 0.0 with 
U = 1.0: p a \ = pal = Pa (solid curves) and p^ = pb± = p a (dashed curves). 

FIG. 3. The J dependence of the band-narrowing factor q in the P-af state. 

FIG. 4. The J dependence of the occupancies in the P-af state: o?2| = c^j. = ^2, t a ^ = t a ^ = t a 
and 4| = tbi = t b . 

FIG. 5. The J dependence of the total energies E of the AF-p state in the GA (bold, solid 
curve) and HFA (thin, solid curve), and E of the P-af state in the GA (bold, dashed curve) and 
HFA (thin, dashed curve). 
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This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0005271vl 



